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Abstract 

The isospin dependence of incompressibility is investigated in the Skyrme Hartree-Fock (SHF) 
and relativistic mean held (RMF) models. The correlations between the nuclear matter incom- 
pressibility and the isospin dependent term of the finite nucleus incompressibility is elucidated 
by using the Thomas-Fermi approximation. The Coulomb term is also studied by using various 
different Skyrme Hamiltonians and RMF Lagrangians. The symmetry energy coefficient of incom- 
pressibility is extracted to be K T = —(500 ±50) MeV from the recent experimental data of isoscalar 
giant monopole resonances (ISGMR) in Sn isotopes. Microscopic HF+random phase approxima- 
tion (RPA) calculations are also performed with Skyrme interactions for 208 Pb and Sn isotopes to 
study the strength distributions of ISGMR. . 

PACS numbers: PACS:21.60.-n, 21.65.+f 
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I. INTRODUCTION 



The study of the correct value of the nuclear matter incompressibility Koo continues to be 
an active area both theoretically and experimentally The value of is the most fundamen- 
tal quantity used in determining the nuclear matter equation of states (EOS). The important 
experimental information on is provided by the isoscalar giant monopole resonances (IS- 
GMR) in finite nuclei. Non-relativistic and relativistic mean-field models have been success- 
fully applied descriptions of many properties of the nuclear ground states and also collective 
excitations including giant resonances. The Skyrme Hartree— Fock (SHF) model is one com- 
monly used non-relativistic mean field model l|, [2j. The Hartree- Fock(HF)+random phase 
approximation (RPA) calculations with Skyrme interactions were also performed to obtain 



the response functions of ISGMR 



SB 



4|. The relativistic mean field (RMF) model is based 



on an effective Lagrangian for the interacting many-body system [5|. The time-dependent 



Hartree (TDRMF) and RPA calculations were performed for the ISGMR {(3, 7] based on the 
same RMF Lagrangian. 

The nuclear matter incompressibility is determined by the second derivative of the 
energy per particle E/A with respect to the density p at the saturation point, 



2 d (h 



(i) 

P — prim 



where h is the isoscalar part of the Hamiltonian density H nm for nuclear matter. The nuclear 
matter incompressibility is not a directly measurable quantity. Instead, the energy of 
ISGMR E isgmr is expressed in terms of the finite nucleus incompressibility Ka as {3} 



/ 1,1 Kx (2) 

L '^"<- W m<r a >J (2) 

where m is the nucleon mass and < r 2 > m is the mean square mass radius of the ground state. 
The finite nucleus incompressibility can be parameterized by means of a similar expansion 
to the liquid drop mass formula with the volume, surface, symmetry and Coulomb terms; 

K A = K OB + K surf A~ 1/3 + K T 5 2 + Kcoui^, (3) 

where 8 = (N - Z) /A. 

In this work, we study correlations among nuclear matter and finite nucleus incompress- 
ibilities for a large number of different Skyrme and RMF parameter sets. In particular, we 



study the correlations between and K T} and also between and Kcoui- Then, we 
extract the values of Koo and K T from experimental ISGMR recently observed in RCNP, 
Osaka University, and also in Texas A&M University. Skyrme HF+RPA calculations are 
also performed to study the detailed structure of ISGMR taking into account the coupling 
to the continuum. In section II, we describe the nuclear matter and finite nucleus incom- 
pressibilities by using the Thomas-Fermi approximation. HF+RPA results of ISGMR will 
be shown for 208 Pb and Sn isotopes in Section III. The summary is given in section IV. 



II. THOMAS-FERMI APPROXIMATION FOR SKYRME AND RMF HAMILTO- 
NIAN DENSITY 



The Skyrme interactions used in many recent studies have nine parameters 
(to, ti, t2, £3, xq, Xi, Xii £3, a) in addition to the parameters of the spin-orbit interaction. The 
energy density functionals of the Skyrme interaction are expressed by using the Thomas- 
Fermi approximation for the kinetic energy density. The nuclear matter properties are 
defined by using these energy density functionals. Various correlations among nuclear 
matter properties have been discussed in the cases of the SHF and RMF models recently 



11] . It was shown that the Skyrme parameters can be expressed analytically in 



terms of the isoscalar and the isovector nuclear matter properties of the Hamiltonian density 
in ref. js]. 

In this section, we detail study of correlations between the nuclear matter incompressibil- 
ity .fToo and the symmetry term K T of the finite nucleus incompressibility, and also between 
and the Coulomb term Kc ou i i n the SHF model for 14 different parameter sets (SI, SIII, 
SJY^Sy^,^y^SkM^k^*^Ly4, MSkA, SkI3, SkI4, SkX, SGI, SGII) taken from Refs. 
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15 



16 



13, 
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19 1 and the RMF model for seven different 



parameter sets 
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22 



23 



24] 



(NL3, NLSH, NLC, TM1, TM2, DD-ME1, DD-ME2) taken from refs. [5|,|20|. 
The non-linear potential of a mesons are introduced in the first five parameter sets of RMF. 
The non-linear potential of the uj meson is also added in the parameter sets TM1 and TM2. 
The meson— nucleon couplings depend on the total density in DD-ME1 and DD-ME2 param- 
eter sets, while the couplings are constants in the other five parameter sets. The non-linear 
potential of a mesons are fixed to be zero in DD-ME1 and DD-ME2 parameter sets. 

The isoscalar part h(p) and the isovector part £<s(p) of the Hamiltonian density H nm are 
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defined by 



h{p) 



lim H n 
lim 



d 2 



H„ 



P 



(4) 
(5) 



2 J->0 dP 

where I is the asymmetry parameter I = (p n — Pp) / P- The derivative terms and the Coulomb 
term in the Hamiltonian density do not give any contribution in the infinite nuclear matter 
calculations. The Thomas-Fermi approximation can be applied for the kinetic energy of 
the Hamiltonian density H nm (p n , p p ) of nuclear matter in the SHF and RMF model. The 
explicit forms of the Hamiltonian densities H nm are found in ref. js]. 

The physical properties of infinite symmetric nuclear matter with p n = p p can be obtained 
from the following six equations: 



o= l( h - 
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J 
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where p nm , E Q , and J are the nuclear saturation density, the binding energy per nucleon, 
the incompressibility of symmetric nuclear matter and the symmetry energy, respectively. 

We denote K T the symmetry term of the finite nucleus incompressibility Ka because the 
symbol K sym has been already used as one of the isovector nuclear matter properties defined 
by Eq. tfTTT) . The volume term of the finite nucleus incompressibility Ka is identified as the 
nuclear matter incompressibility K^. The symmetry contribution K T is related to nuclear 



matter properties as 



3], 



K T — K S y m + 3L 



27Lp 2 nm d*h 
Koo dp 3 
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The analytic formulas for K sur f and Kcoui are given by 



K SU rf = 4tt r 2 



4cr(p nm ) + 9p 

ft? 



dp 2 



+ 



hAa(p nm )p 2 nm d 3 h 
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dp 3 



P=Pn 
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where r is the radius constant defined by 



r 



Anp ri 



1/3 



(15) 



In Eq. (1131) . a is a surface tension in symmetric semi- infinite nuclear matter defined by 

h(p, 



H(p)~ 



Pr. 



-P 



dz, 



(16) 



where H is the Hamiltonian density. K sur f can be evaluated by the extended Thomas- Fermi 
approximation and the scaled HF calculations on semi-infinite nuclear matter in the SHF 
model. These evaluations show an approximate relation K sur f ~ — Koo within an accuracy of 
a few % in the SHF model. In RMF, the study of an extended Thomas-Fermi approximation 
gives a slightly larger surface contribution, for example, K sur f ~ —1.16^00 in the case of 
NL3. 

The values of K T and Kcoui are calculated by using various Skyrme Hamiltonians and 
RMF Lagrangians and shown in Figs. [HandEJ K T is largely negative and has anti-correlation 
with the nuclear matter incompressibility K^. Namely, any Hamiltonian which has a larger 
gives a smaller K T . The variations of K T for the Skyrme interactions are 



K T = (-400 ± 100) MeV 



for Skyrme interaction. 



(17) 



On the other hand, the values of RMF are largely negative and have more variation among 
the seven effective Lagrangians, 



K T = (-620 ± 180) MeV 



for RMF Lagrangian. 



(18) 



In principle, the value Kc ou i should be model- independent. Nevertheless, we can see a weak 
correlation between and Kc ou i in Fig. [2J The correlation between and Kc ou i can 
be expressed analytically by using the isoscalar nuclear matter properties as given in the 
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Appendix. Among the 13 parameter sets of Skyrme interactions, the variation of Kcoui is 
rather small, 

Kcoui = (-5.2 ± 0.7) MeV 

compared with that of K T . The values of Kcoui m RMF show essentially the same trend, 
but have a larger variation. 

Recently, the ISGMR strength distributions in the Sn isotopes from 112 Sn to 124 Sn have 



been measured by using inelastic a scatterings at RCNP, Osaka University 30j. The ISGMR 



of Sn isotopes were also studied in Texas A & M University 



26 



29J. The value of K T was 



extracted to be K T = (—395 ± 40) MeV for Sn isotopes assuming Ka to be a quadratic 
relation with the symmetry parameter 6 = JV ~ Z , i.e., Ka=c+K t 5 2 where c is a constant. 
We should be careful to apply this relation to obtain the value K T since the surface and the 
Coulomb contributions in Eq. ([3]) are also functions of the mass number A. We examine 
the K T dependence of the E ISGMR by using the formula ([3]). In Fig. 3, the difference of the 
compressibility AKa = K A — Ka=h2 is plotted as a function of 5 = (N — Z) / A. We adopt 
four Skyrme interactions and two RMF Lagrangians. In the analysis, the surface term is 
taken to be K SUJ .f = —K^ for Skyrme model and K sur f = —lASK^ for RMF model. The 
adopted interactions vary from a smaller K T value of -350 MeV for SkM* to a larger value 
of -700 MeV for NL3. The empirical isospin dependence of AKa is close to the results of 
SIII, SIV and DD-ME1, which have K T = -(500 ± 50) MeV. The present extracted value is 
about 30% larger than the value reported in ref. [30j|. This difference is mainly due to the 



mass number dependence of K sur f and Kcoui which were neglected in ref. [301 ] . 



III. HF+RPA CALCULATIONS FOR ISGMR 

In order to extract the incompressibility for infinite nuclear matter from the exper- 
imental values of the ISGMR, self-consistent RPA calculations were performed and results 
thereof were compared directly with the experimental data both in SHF and RMF model 
. The sum rule approach with the constrained HF calculations was also adopted 



to disentangle in ref. 25j. The experimental data of 208 Pb was adopted in these studies 
since ISGMR data is most well-established in this nucleus. There has been an attempt to 
extract the volume, surface, symmetry, and Coulomb terms in Eq. ([3]) by using various sets 
of experimental data of ISGMR in different nuclei j^]. However, the results depended very 
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much on the adopted data set and gave a negative conclusion to the feasibility of simulta- 
neously determining the four terms of Eq. ([3]) from experimental data. In this study, we 
perform HF+RPA calculations for 208 Pb, 112 Sn, 116 Sn, 120 Sn and 124 Sn to enable quantitative 
comparisons to be made with the experimental strength distributions of ISGMR and try to 
extract realistic values of and K T . 
The RPA strength distribution 

S(E) = I < n\Q\0 > \ 2 6(E - E n ) (19) 

n 

is calculated by using the IS monopole operator 

Q A=0 ' r=0 = ^E^ ( 2 °) 



'Ait 

■i 

The k-th energy moment of the transition strength is defined by 

m k = J dEE k S(E). (21) 

The average energy will be obtained by a ratio between the moments mi and m , 

E = m 1 /m . (22) 

The scaling model of ISGMR gives the excitation energy 

E s = >/m 3 /mi, (23) 

while the excitation energy produced by the constrained HF model is written as 

E c = y/mjm^. (24) 

The excitation energies defined by Eqs. (!22|) - (|24|) are identical in the case of a sharp single 
peak which exhausts 100% of the sum rule. However, in reality, both the experimental data 
and the calculated results show a large width of a few MeV even in the most well-established 
ISGMR in 208 Pb. In particular, the scaling energy E s has a large uncertainty due to the 
high energy tail of monopole strength, which is always the case in experimental data. On 
the other hand, E and E c are rather close within a 0.1~0.2 MeV difference even when the 
ISGMR peak has a large width. Because of these reasons and also due to the theoretical 
clear background, we identify E c as Ejsgmr to obtain Ka through Eq. (2). It should be 



noticed that the rms mass radius is needed to obtain Ka through Eq. (2), but not the charge 
radius. Since the mass radii of Sn isotopes were not determined experimentally so far, we 
use HF mass radii to analyse Ka- 

The continuum HF+RPA results for ISGMR in 208 Pb are shown in Fig. 4 with three 
different Skyrme interactions SIII, SGI and SkM*. The incompressibility of The SIII inter- 
action is 355 MeV, while those of SGI and SkM* are 256 MeV and 217 MeV, respectively. 
Various excitation energies ff22l . (T23]) and fT24"|) are given in Table [H As far as the excitation 
energy is concerned, SkM* shows a fairly good agreement with the experimental data ob- 
tained by (a, a') scatterings [26]. The average energies E of SGI and SIII are 1.0 MeV and 
3.3 MeV higher than the empirical one, respectively. The calculated width for the SkM* in- 
teraction shows almost the same width as that of experimental data. This agreement implies 
that the dominant contribution of the width of ISGMR stems from the Landau damping 
and the coupling to the continuum, which are properly taken into account in the present 
calculations. On the other hand, the coupling to the many-particle many-hole states might 
have a minor effect on the width of ISGMR. The agreement between the results of SkM* 
and the experimental data strongly suggests a lower incompressibility Koo ~220 MeV as a 
realistic one for infinite nuclear matter. The two-body spin-orbit and the two-body Coulomb 
interactions are not taken into account in the present RPA calculations although the HF 
calculations include both interactions. It was pointed out in ref. 3JJ that the net effect 
of the two interactions in RPA decrease the centroid energy of ISGMR E ISGMR in 208 Pb 
by about 300 keV. Then, an approximate relation (SK^/ Koq) =2(5Ejsgmr/^isgmr) from 
Eq. (2) gives a decrease of the incompressibility by dKoo ~10 MeV. Due to this effect, a 
realistic incompressibility should be slightly larger than than of SkM* to be ~230 MeV. 



In ref. 



251 ]. the nuclear matter incompressibility Koo was discussed taking into account the 



full interactions of various Skyrme parameter sets in the constrained HF model. Comparing 
the calculated excitation energies E c in Eq. (T241) with the experimental data of ISGMR in 
208 Pb, they concluded that the nuclear matter incompressibility should be i^oo=230-240 MeV 
in order to explain the experimental data. This conclusion is consistent with the present 
HF+RPA calculations with SkM* since the two-body spin-orbit and Coulomb interactions 
give the net effect of SK^ ~10 MeV. 

Next, let us discuss the isospin dependence of the excitation energies of ISGMR in Sn 
isotopes. For the HF+RPA calculations of Sn isotopes, we adopt the SkM* interaction, 
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TABLE I: The energies of ISGMR in 208 Pb and the finite nucleus incompressibility Ka calculated 
by using Eq. ([3]) with Ejsgmr = E c . The HF+RPA results are calculated by using SkM*, SGI 
and SIII interactions. The sum rules are obtained by summing up the strength up to E x .=25.5 
MeV. The average energy is defined by the ratio E = mx/rriQ in Eq. (|22p . while E s and E c are 
obtained as E s = \/m^/mi and E c = \/mi/m-\ in Eqs. (|23p and (|24p . The experimental data 



are taken from ref. 



261 ] . The experimental value Ka is calculated by using E c and the mass radius 



of the SkM* interaction. All values are given in units of MeV. 





E 


E s 


E c 


K A 


#00 


SkM* 


13.87 


14.29 


13.79 


141.5 


217 


SGI 


15.17 


15.47 


15.07 


170.0 


256 


SIII 


17.47 


17.88 


17.31 


226.4 


355 


exp (Texas A&M) 


14.17±0.28 




14.18±0.11 


149.6±2.3 





discretizing the continuum with a large harmonic oscillator basis up to the maximum major 
quantum number N max =lQ. The RPA calculations are performed by using the filling ap- 
proximation for the neutron orbits, i.e., the neutrons occupy the orbits from the bottom of 
the potential to the Fermi level in order. The last neutron orbit is partially filled according 
to the neutron number. The calculated values of the transition strength B(E) are averaged 
by using a weighting factor p{E x — E) 



S(E x )=Y,B(E i )p(E x -E l 



(25) 



where 



p{E x - Ei 



A/2 



(26) 



7T (E x - Ei) 2 + (A/2) 2 ' 
The width parameter is taken to be A=l MeV in Fig. [5j We also performed the continuum 
RPA calculation in a nucleus 116 Sn and found essentially identical results as far as the 
excitation energies listed in Table [Til are concerned. As seen in Fig. [5j the RPA results 



show reasonable agreement with the experimental data obtained by Texas A&M [26|, |29] 
and also by RCNP [30|. The various average energies E = mi/m 0} E s = ^Jm^/mi and 
E c = ypm\[mZ\ are listed in Table [UJ In general, the calculated average energies by RPA 
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with SkM* are few hundreds keV higher than the empirical ones. The two sets of recent 
experimental data show slight differences in 112 Sn and 124 Sn which should be confirmed in 
the future by further experimental study. The calculated excitation energies decrease from 
112 Sn and 124 Sn by about 1 MeV which is consistent with the observed data. This decrease 
is expected from a large negative symmetry term K T in the finite nucleus incompressibility 
discussed in Section 2. 

The RPA results give about 80% of the observed widths of ISGMR in Sn isotopes. These 
results are almost identical to the case of 208 Pb. Thus, we can conclude from the results of 
Sn isotopes and 208 Pb that the major part of the width of ISGMR stems from the Landau 
damping and the coupling to the continuum. 



TABLE II: The energies of ISGMR in Sn isotopes. The HF+RPA results are obtained by using 
the SkM* interaction. The RPA sum rules are obtained by summing up the strength up to 
E x =25.5 MeV to be consistent with experimental data. The average energy obtained by the ratio 
E = rrii/rriQ in Eq. (|22p . while E s and E r ar e defined by E s = \f / m~ and E c = \Jm\jm^.\. The 
experimental data are taken from refs. [2a, |29J (Texas A&M) and 30](RCNP). The experimental 



29| (Texas A&M) and 

value Ka is calculated by using E c and the mass radius of the SkM* interaction. All values are 
given in units of MeV. 





RPA(SkM*) 


Texas A&M 


RCNP 




E 


E s 


E c 


K A 


E 


E s 


E c 


E 


E s 


E c 


K A 


112 Sn 


17.1 


17.5 


16.9 


142.7 


15.43 


16.05 


15.23 


16.2 


16.7 


16.1 


129.5 ±1.6 


116 Sn 


16.6 


17.0 


16.5 


140.2 


16.07 




15.90 


15.8 


16.3 


15.7 


126.9 ±1.6 


120 Sn 


16.1 


16.5 


15.9 


133.7 








15.7 


16.2 


15.5 


127.0 ±1.7 


124 Sn 


16.2 


16.7 


16.1 


137.1 


14.50 


14.96 


14.33 


15.3 


15.8 


15.1 


122.7 ±1.7 



IV. SUMMARY 

In summary, we studied the ISGMR of 208 Pb and Sn isotopes in order to disentangle the 
nuclear matter incompressibility and also the symmetry term K T in the finite nucleus 
incompressibility. Firstly, the Thomas Fermi approximation is adopted to obtain the various 
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terms of the finite nucleus incompressibility Ka in the SHF and RMF models. The corre- 
lations between and K T and also between and Kcoui are elucidated in the various 
sets of the Skyrme Hamiltonians and RMF Lagrangians. We extracted the symmetry term 
to be K T = —(500 ± 50) MeV from the analysis of the isospin dependence of the excitation 
energies of ISGMR in Sn isotopes. Secondly, we perform HF+RPA calculations to study 
detailed structure of ISGMR in 208 Pb and Sn isotopes. In 208 Pb, the results of continuum 
RPA with three sets of Skyrme interactions SkM*, SGI and SIII are compared with the 
experimental data. It is shown that SkM* gives a satisfactory description of ISGMR of 
208 Pb in both the excitation energy and the width. The nuclear matter incompressibility 
is extracted to be ~230 MeV from the RPA analysis with the additional effects of the 
two-body spin-orbit and Coulomb interactions. The isospin dependence of ISGMR in Sn 
isotopes is also studied by using the HF+RPA calculations with the SkM* interaction. We 
pointed out that the RPA results give a reasonable account of the isospin dependence of the 
observed data, although the calculated excitation energies are few hundreds of keV higher 
than the observed ones. It should be noted that appreciable differences exist between the 
two measurements in the excitation energies of 112 Sn and 124 Sn, listed in Table [TTJ, which 
makes it difficult to accurately determine the value K T of the finite nucleus incompressibility. 
This problem remains as a future challenge both in experimental and theoretical studies. 
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TABLE III: The values of a and b in Eq. (A. 4) as a function of the effective mass m*. 



Tfl* J TYl 


Qj 


b 


0.6 


-0.0104314 


-2.15499 


0.7 


-0.0105882 


-2.10745 


0.8 


-0.0106691 


-2.08292 


0.9 


-0.0107185 


-2.06796 



APPENDIX: ANALYTIC FORMULAS FOR K Coui 



The Coulomb term Kq uI can be expressed by nuclear matter values and p nm , and 
the Hamiltonian density for symmetric nuclear matter h as follows 



Kcoul = 3 ^ ( 1 27 Pnm d " h 

( '"" 5r I dp 3 



P — pnm 



5 r V K, 



27 



■— cp^(3a + 1) - -(a + 1)^0 + 



(A.l) 



where c = (37r 2 /2) 2 ^ 3 and a is the power of the nuclear density-dependent term in the 
Skyrme interaction |l|]. The nuclear mater incompressibility Koo is also expressed as 



10m 



cp 



2/3 



3(3a- 1) -2(3a-2) 



m 
m 



+ 9(a + l)E . 



(A.2) 



The linear power of the density-dependent term of the Skyrme interaction a can be elimi 
nated from Eq. flA.lj) by using Eq. flA.21) ; 

3e 2 



K, 



Coul 



5r 

x 



^ + 3r? ( -1 + \m-}/) - 2E + J5- |^ (1 - m.//) + E (27 - 25™^) 



^ (3 - 2m-} f ) + E _ 

2^ 2/3 



5 
-1 



(A.3) 



where m eff = m*/m and r) = h 2 cp^/m. Eq. flAj) is proportional to the power of 1,0 and 
-1 of i^oo. However, this equation can be well parameterized to be 



K, 



Coul 



aK^ + b, 



(A.4) 



where a and b depend on the effective mass m*. With the standard value of the saturation 
density and the energy of nuclear matter p nm =0.16 fm~ 3 and E =1Q MeV, the values of a 
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and b are tabulated in Table Al. 
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FIG. 1: Correlation between the nuclear matter incompressibility Koq and the symmetry contri- 
bution K T in the finite nucleus incompressibility calculated by using 14 parameter sets of the SHF 
(open circles) and seven parameter sets of RMF (filled circles). The numbers denote the different 
parameter sets: 1 for SI, 2 for SIII, 3 for SIV, 4 for SVI, 5 for Skya, 6 for SkM, 7 for SkM*, 8 for 
SLy4, 9 for MSkA, 10 for SkI3, 11 for SkM, 12 for SkX, 13 for SGII, 14 for SGI, 15 for NLSH, 16 
for NL3, 17 for NLC, 18 for TM1, 19 for TM2, 20 for DD-ME1 and 21 for DD-ME2. 
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FIG. 2: Correlation between the nuclear matter incompressibility and the Coulomb contribu- 
tion Kcoul in the finite nucleus incompressibility calculated by using 14 parameter sets of the SHF 
(open circles) and seven parameter sets of RMF (filled circles). See the caption of Fig Q] and the 
text for details. 
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Experimental data are determined by using the excitation energies of ISGMR in ref. 
HF mass radii. See the text and the caption to Table II for details. 
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FIG. 4: (Color online) Continuum HF+RPA response functions of Skyrme interactions with the 
different incompressibilities SIII, SGI and SkM*. The experimental data are taken from ref. 
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FIG. 5: (Color online) HF+RPA response functions of 112 Sn, 116 Sn, 120 Sn and 124 Sn nuclei with 



Skyrme interaction SkM*. The experimental data are taken from refs. 
(RCNP). 
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